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Abstract. The developed scientific works are directly related to solving the problems of identifying 
multi-parameter objects. The proposed algorithm is effectively implemented on a computer, which is 
confirmed by calculations of problems of supersonic aerodynamics. 
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I. INTRODUCTION 


Let us consider the methods of identification and assessment of the state in solving problems of 
the process of flowing around an oscillating profile by a supersonic flow. 

Using the proposed algorithm, we obtain identification models. As an aerodynamic characteristic 
used shock wave departure A . The main factors influencing the response function are: X_1-Mach number 
M_ o, X 2-relative distance from the bottom edge of the profile to the screen H , X_3-angle of attack of 
the profile a, X_4-relative thickness of the profile C. The process of building the model consists of three 
stages. At the first stage, the linear part of the model is considered. The second is the quadratic part and 
the third is the nonlinear part. 

To achieve the resulting conditioning of the matrices, the experiments were conducted by the 
method of rotationally central compositional planning. 

At the first stage, the linear model has the form 


rl a 

ya + di-1 GX; 
On the first row of selection we get six combinations (of four in two) linear polynomials. The private 
polynomial can be written as: 


43; = Aoi + a,;X1 + Az; Xj44,1 — 13 


V2 i42 = Ap 43 + Ay j43X2 + Az j43Aj42,1 = 12 





Yo4 — og + A16X3 + AzgX a4. 


To determine the coefficients a_0i, a_li, a_2i, (i = (1,6) ) we compose a system of conditional equations 
in the form 


VN wk vN Ny _ 
k=14 ij Xi — Qo deena Xix — Gi digai Xie — G Liar XjeXix = 0 
k=1 ijAjk — A dient Aje — Ui Lika XinXjx — G dia A jk =0 


i=1,n,j=itin 


(1) 


where 
The system of equations is solved by the Gauss method. The standard error is calculated by the formula 


s2 lyn vy Po\2 ._ 4 2 
Os ~~ 5 uk=1( ‘ik .) = 5, Cr 
ye | ee | , a 
where ~ J& — experimental values at the i-th point; ~ “/* — estimated value of private descriptions; 


(2) 


2 
O5 — standard error of the sequence of models. 


Il. MAIN PART 


We will gradually increase the complexity of the partial and complete regression equations 
and at the same time trace the nature of the change in the mean square error. 
The first stage of building the model — on the first row of selection, the following coefficients 
were obtained, and accordingly, the standard deviations (Table 1). 
Table. 1. 
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52, = 0.4728 


In the next series of selection, those variables that satisfy the condition 


co? oC ¢ : —— 
6; < bm i= 1,6 

Among the obtained variables satisfying the criterion U_14, U_23, U_34. 

In the second row of selection, we obtained with a combination (of three by two) linear polynomials 


214-23 = Boy + Dy Vi4 + b21¥23, 244-34 = Boz + by2V14 + bo2¥34, 
£93-34 — bo3 + by3¥23 + b23Y34. 
The numerical values of the table 2 are obtained. 
Table. 





bf = 0.3298 


From the obtained results it is seen that only Z_ (14-23) satisfy criterion (2). 
Then Z14-23 = Do + 04414 + 91923 
Via = A93 + Ay3Xy + Ag3X4 Yo3 = Apg + AygXq + Az4Xz 
244-23 = boy + by1A93 + b21Ao4 + by 4A43Xq + B21 A44Xq + Dg1Az4X3 + by A23Xq 
Ag = bo, + b44Q03 + b214Ao94 = 1.4714, a, = b,,a,3 = —0.4614 
dy = by1d44 = 0.0696, a3 = by dy, = —0.5224, ay = 0.063 


At the first stage, a linear model is obtained. 


y= 1.471 — 0.4614X, + 0.0696X, — 0.5224X3 + 0.063X,_ 


The second stage of building the model. 
In the first row of selection, the following coefficients were obtained, respectively, with standard 
deviations (Table 3). 


Table 3. 





56) = 0.2479 


In the second row of selection (table. 4.) 


a ee 0.9265 0.2423 





56 = 0.2425 
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In the third row of selection 
iq — 0, Ci, — —0.1206, Cyo — 1.1182, 6* = 0.2418 
V1214-1324 = Cor + 1212-14 F C21 213-14, 
212-14 — 0924012912 z b59¥14,213-14 = bo3 * bi3¥13 + b53¥ 14 
Viag 6 $623 6 KH 1 


Vi214-1314 = C11D12A91 + Cy1b22A93 + C21b13402 + (G1 Di2Qy1 + C41 22013 + 
€71b13013 + C21b23013)Xq + (Cy1by2491 + C11 b22A93)X2 + Cy1b1323X3 + 
€21b93073X4 


At the second stage, the model is obtained 
Y = 0.0224 — 0.0617X2 + 0.014X5 — 0.0002X2, 62 = 0.2418. 
The third stage of building the model. 
In the first row of selection we obtain the following coefficients, respectively, with standard deviations 
(Table 5). 


Table. 5. 





dg = 0.2310 


In the second row of selection 


bp, = 0.0921, b,, = —6.1055, 67 = 0.195 


At the third stage, the model is obtained 

Y = 0.0558 + 0.0213X,X3 + 0.0085X3X, — 0.005X,X3X, 

And so the object identification model 

Y = 1.549 — 0.4614X, + 0.069X, — 0.5224X3 + 0.063X, — 0.0617X? + 0.014X2 — 
0.0002X2 + 0.0213X>X3 + 0.0085X3X, — 0.005X,X2X,. 
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Fig. 1. Root mean square error from the 
number of breeding rows 
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Fig. 2. Root mean square error on the number of process additivity 


The dependence of the mean square error on the number of selection series in the first (curve 1), second 
(curve 2), and third stages (curve 3) are shown in Fig. 1, and the dependence of the mean square error on 
the additivity number is shown in Fig. 2. Thus, a methodology was developed for constructing nonlinear 
multipara meter mathematical models and was used in modeling and identifying parameters of the 
process of a supersonic flow around an oscillating profile. 
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